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ABSTRACT 

Recently, Widrow and collaborators announced the discovery of vertical density 
waves in the Milky Way disk. Here we investigate a scenario where these waves were 
induced by the Sagittarius dwarf galaxy as it plunged through the Galaxy Using nu- 
merical simulations, we find that the Sagittarius impact produces North-South asym- 
metries and vertical wave-like behavior that qualitatively agrees with what is observed. 
The extent to which vertical modes can radially penetrate into the disc, as well as their 
amplitudes, depend on the mass of the perturbing satellite. We show that the mean 
height of the disc is expected to vary more rapidly in the radial than in the azimuthal 
direction. If the observed vertical density asymmetry is indeed caused by vertical os- 
cillations, we predict radial and azimuthal variations of the mean vertical velocity, 
correlating with the spatial structure. These variations can have amplitudes as large 
as 8 km . 

Key words: Galaxy: disc, structure - galaxies: formation - galaxies: kinematics and 
dynamics - methods: A^-body simulations 



1 INTRODUCTION 

Minor mergers can significantly perturb the overall struc- 
ture of their host galactic disc ( Quinn et al.|p!993 Villalo- 
|bos &; Helmi|2008J . As they merge, relatively small satellite 
galaxies can induce the formation of spiral arms and ring- 
like structures as well as radial migration, significantly fiare 
or warp the disc, and infiuence the growth of a central bar 
([Kazantzidis et al. 2008; Younger et al. 2008 ; Minche v et al. 
2009,.Quillen et al. 2009, Purcell et al.^2011; Bird et al.|2012| 



by this interaction. An example is the kinematically cold 



Gomez et al.|2012b 



Purcell et al. 



(2011 



hereafter Pll) presented simula- 
tions of the response of the Milky Way (MW) disc to tidal 
interaction with the Sagittarius dwarf galaxy (Sgr). They 
showed that many of the global morphological features ob- 
served in the Galactic disc can be simultaneously explained 
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structure known as the Monoceros ring ( Newberg et al. 



2002 



p'uric et al.^^2008 ), which naturally emerges in these simula- 



tions, although its origin is still a matter of debate (see [Conn 
et al.||2012|. I Gomez et al.l (|2012a" hereafter G12) showed 



that perturbations observed in the phase-space distribution 
of old disc stars in the Solar Neighbourhood (SN) can be 
reproduced qualitatively with these simulations. They in- 
terpreted such perturbations as signatures of radial density 
waves excited on the plane of the disc by Sgr. Perturbations 
in the vertical direction were not explored in their work. 
However, using a much larger photometric and spectroscopic 



data set, Widrow et al. (2012 hereafter W12) recently iden- 
tified a North-South asymmetry in both the spatial density 
and the velocity distribution of SN stars. The asymmetry 
has the appearance of a coherent, wave- like perturbation, 
intrinsic to the disc. W12 speculate that this perturbation 
could have been excited by the passage of a satellite galaxy 
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through the Galactic disc. In this Letter, we explore the pos- 
sibility of Sgr being the perturber associated with both the 
vertical and radial modes. 



2 SIMULATIONS 

In this Section we briefly describe our simulations; we re- 
fer the reader to Pll and G12 for a more detailed descrip- 
tion. Two simulations with different models for the Sagit- 
tarius dwarf galaxy progenitor were performed. A Light 
(Heavy) Sgr progenitor, with effective virial mass Mvir = 
10^°-^ M© (10^^ Mo), was initialised with a NFW dark mat- 
ter halo of scale length 4.9 kpc (6.5 kpc), self-consistently 
with a separate stellar component. For the stellar compo- 
nent, a King profile with core radius 1.5 kpc, tidal radius 
4 kpc, and a central velocity dispersion of 23 km s~^ (30 
km s~^) was used. The satellites were launched 80 kpc from 
the Galactic centre in the plane of the MW, traveling ver- 
tically at 80 km s~^ toward the North Galactic Pole. The 
mass loss that would have occurred between virial radius 
infall and this "initial" location is accounted for by trun- 
cating the progenitor DM halo mass profile at the instanta- 
neous Jacobi tidal radius, n = 23.2 (30.6) kpc. This leaves 
a total bound mass that is a factor of ^ 3 smaller than 
the effective virial mass originally assigned. The simula- 
tions reach a present-day configuration after ^ 2.7 Gyr (2.1 
Gyr) of evolution. In both cases, the host galaxy includes 
a NFW dark matter halo with a scale radius Ts = 14.4 kpc 
and virial mass Mvir = 10^^ M©; the disc has a mass of 
3.59 X 10^° M©, an exponential scale length of 2.84 kpc, 
and a vertical scale height of 0.43 kpc; the central bulge has 
a mass of 9.52 x 10^ M© and an n = 1.28 Sersic pro- 
file, with an effective radius of 0.56 kpc. The initial disc in 
both Sgr-infall models is completely smooth at t = Gyr. 
The simulations followed the evolution of 30 million particles 
with masses in the range 1.1 - 1.9 x lO^M©. 
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Figure 1. Overdensity map of the Heavy (left) and Light (right) 
*S'5'r simulations at the present time. The colours and relief indicate 
the ratio of the local stellar density to the mean axysimmetric disc 
density at the corresponding galactocentric distance. The blue 
and red dots indicate the location of the volumes explored. 
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Figure 2. Left panel: Total (dashed lines) and in-plane (solid 
lines) energy distributions of disc particles located within the blue 
and red volumes shown in Fig^ Right panel: As in the left panel, 
but for the vertical energy distribution. 



3 PERTURBATIONS IN LOCAL VOLUMES 

Fig. [l] shows an overdensity map of the Heavy (left) and 
Light (right) Sgr simulation discs at present-day configura- 
tion. The maps are obtained by normalising the local stel- 
lar density to the mean axysimmetric density at the cor- 
responding galactocentric distance. The non- axysimmetric 
energy kick imparted by the satellite as it merges with the 
host induces the formation of spiral density waves. As shown 
by G12, waves excited by the Heavy Sgr satellite can be 
detected in SN-like volumes as peaks in the local energy 
distribution. These density waves are mainly in the radial 
direction. In fact, in very small local volumes (i.e., distances 
^ 0.2 kpc), these waves can be observed as well-defined fea- 
tures in the radial (vr) and tangential (i;^) velocity field (see 
Minchev et al.|2009| G12). In the left panel of Fig.|2] we show 
with dashed lines the normalised total energy distribution 
E = /(x, v) of two SN-like cylindrical volumes extracted 
from the Heavy Sgr simulation. The volumes have a 1 kpc 
radius, and are located at 8 kpc from the galactic centre. 
Their locations, chosen based on G12 results, are indicated 
with a blue and a red dot in the left panel of Fig.^ Note the 
well-defined peaks in these distributions, which reveal the 



presence of density waves. As expected from disc orbits, the 
in-plane energy distribution E\\ = /(x, t'r , t'c/,) (solid lines) 
is very similar to the total energy distribution, indicating 
that these peaks are associated with perturbations mainly 
in the plane of the disc. Nevertheless, the stellar particles in 
these volumes present non-symmetrical energy distributions 
in the vertical direction Ez = /(x, Vz), as shown in the right 
panel of Fig|2] Furthermore, these distributions are shifted 
with respect to one another. 

To explore whether the asymmetries and shifts observed 
in the Ez distributions are indications of vertical density 
waves, we compute for both volumes the distribution of stel- 
lar particles, as a function of height with respect to the 
midplane of the galactic disc, n(Z)|^^, or simply n(Z). A 
north-south asymmetries in this distribution could be an 
indication of vertical modes (see W12). To compute n(Z), 
we have carefully aligned the disc with the X-Y plane. This 
is done by iteratively computing, and aligning with the Z- 
direction, the total angular momentum of the disc particles 
located within 4 kpc radius cylinders of decreasing height. 
In order to identify signatures of vertical density waves, it is 
desirable to compare n(Z) obtained from the perturbed disc 
with its corresponding smooth underlying distribution. For 
this purpose, W12 fitted a smooth two-component model to 
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Figure 3. First row: The blue and red dots show the distribution 
of particles as a function of height, n(Z), with respect to the 
midplane of the disc, obtained from the volumes shown in Fig^ 
The black solid lines show a model of the underlying smooth 
distribution nav(Z). Second row: Distribution of the residuals, 
A = {n{Z) — nav(Z))/nav(Z), for both volumes. Third row: Mean 
vertical velocity as a function Z. Error bars indicate Poisson noise. 



their stellar sample's number density. In this work we fol- 
low a different approach: we obtain a smooth distribution 
of stellar particles, as a function of height, by azimuthally 
averaging n(Z) i.e., 

nav(Z)|^ = (2^)-i r n{Z,0)^^ d9. 



Our assumption is that local asymmetries of this distribu- 
tion are erased after averaging over all azimuthal angles. In 
addition, we expect nav(Z)|^ to be a better representation of 
the true smooth height distribution of particles than smooth 
analytic fits. The top panels of Fig. |3] show, with coloured 
dots, the n(Z) distributions obtained from the "blue" and 
"red" volumes in Fig. [l] whereas the black solid line corre- 
sponds to nav(Z). Note that nav(Z) was re- normalised to the 
total number of particle in each volume. Due to finite par- 
ticle numbers, we are able to reliably track n(Z) only up to 
|Z| ^ 1.4 kpc. These panels indicates a shift of the local with 
respect to the smooth distribution, although this occurs in 
different directions for the two cases. Note, however, that 
for large \Z\ (i.e. > 0.7 kpc), the shift in both distributions 
becomes progressively smaller, exhibiting a wave-like pat- 
tern. As in W12, we plot the residual, A = ^^^^ — nav(Z) ^ 

nav(Z) 

to highlight these asymmetries. This is shown in the sec- 
ond row of Fig. |3] In both cases, A is an odd function of 
Z. At least within |Z| < 1.4 kpc, the residuals associated 
with our "blue" volume resemble the distribution presented 
in Fig. 1 of W12, which was obtained from a sample of stars 
in the SN. The amplitude and extent of the perturbations 
found in this simulation are larger than those observed in 



the SN. This indicates that either the Heavy Sgr model is 
too massive or that we are looking at the perturbation at 
an earlier stage. In the bottom row of Fig. [s] we show (-^z) 
as a function of Z. We note that results associated with 
this distribution should be taken with cautious due to their 
relatively low statistical significance. While the "red" vol- 
ume only presents departures from (vz) =0 at large Z, the 
"blue" volume shows a global trend towards negative values 
as Z decreases. As explained by W12, this suggests a coher- 
ent motion of stellar particles away from the disc midplane. 
Note again that the distribution of (vz) associated with the 
blue volume exhibits good agreement with the distribution 
shown in Fig. 4 of W12. 

It is important to realise that the phase-space distri- 
bution of the blue SN-like volume qualitatively reproduces 
not only the observed vertical structure (W12), as we just 
demonstrated, but also features seen in the plane of the lo- 
cal Galactic disk: As shown by G12, local samples of old 
disc stars present a total energy distribution that is in good 
agreement with the distribution shown in the left panel of 
Fig. [2] (blue lines). 



4 AZIMUTHAL AND RADIAL DEPENDANCE 

We have shown that the phase-space distributions of So- 
lar Neighbourhood-like volumes extracted from our Heavy 
Sgr simulation exhibit wave-like perturbations in the ver- 
tical direction. However, it is not clear from the analysis 
shown thus far whether these are localised perturbations or 
are signatures of a global mode perturbing the entire disc. 
From Fig.js] we know that the n(Z) distributions of the two 
analysed volumes are shifted with respect to one another, 
as well as with respect to nav(Z). If the observed shifts are 
signatures of vertical waves, we would expect to find cor- 
relations in the mean height of the disc (Z) as a function 
of galactocentric radius and azimuthal angle. We explore 
this in Fig. [4] where we show maps of (Z) for the Heavy 
and Light Sgr simulation disc out to 20 kpc. To ob- 

tain these maps, we grid the disc with a regular Cartesian 
mesh of bin size = 0.5 kpc aligned with the X-Y plane. On 
each grid node we centre a 1 kpc radius cylinder and com- 
pute (Z) by fitting a Gaussian distribution to n(Z). In the 
top left panel, we show the (Z)-map reflecting the initial 
conditions of both discs at t = Gyr. As expected from 
an initially unperturbed disc, the maps are consistent with 
(Z) = at all radii, except for small-scale departures due to 
finite particle resolution. The situation dramatically changes 
when we explore the present-day configurations. The second 
column of panels in Fig. [4] shows the (Z)-map of the Heavy 
(top) and Light (bottom) Sgr simulation's discs, obtained 
after t ^ 2.1 and 2.6 Gyr, respectively. The colour cod- 
ing indicates different values of (Z) and white regions repre- 
sent volumes devoid of particles. We can now clearly identify 
wave-like, spiral perturbations traveling across the disc. It is 
possible to appreciate how the mean height of the disc grad- 
ually increases (decreases) as we follow one of these patterns 
azimuthally. Departures of (Z) with respect to the midplane 
can be as large as 0.3 kpc, especially for the Heavy Sgr sim- 
ulation (top panels). The local volumes analysed previously 
are indicated with coloured dots. As expected from Fig. [Sj 
the blue (red) volume is located in a region lying slightly 
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Figure 4. First column: Map of the simulated galactic disc's mean height, (Z), at t = Gyr, used in both simulations out to = 20 
kpc. Second column: Maps of the simulated galactic disc's (Z) at present-day configuration, obtained from the Heavy (top) and Light 
(bottom) Sgr simulations, respectively. The different colours and relief indicate different values of (Z) in kpc. Third column: Different 
colours indicate different values of (Z) in kpc, whereas the relief traces the overdensity maps shown if Fig. ^ Fourth column: As before, 
but for the mean vertical velocity (in km s"-*^). The coloured dots in the top panel indicate the location of the volumes analysed in Sec.ls] 



below (above) the midplane. In this simulation, the inner 7 
kpc of the disc has not been strongly vertically perturbed 
by the satellite galaxy. Although the Light Sgr simulation 
(bottom panels) also exhibits wave- like perturbations, they 
generally have smaller amplitudes. Note that strong pertur- 
bations are observed only at > 10 kpc. In both cases, such 
boundary is located near corrotation, where density waves 
are reflected backwards to the outer disc. This can be more 
clearly appreciated in panels A and B of Fig. [S] where we 
show the (Z)-maps in polar coordinates. The black dashed 
lines, at = 8 kpc and = 1.06 rad (measured with re- 
spect to the positive X semi-axis), cross at the location of 
the blue volume. For comparison, this location is indicated 
in both discs. Notice that the inner regions of the Light Sgr 
disc presents a value of (Z) ^ for a larger radial extent. 
Panel C in Fig. [s] shows (Z), as a function of ^, at = 8 
kpc for both discs. Significant vertical perturbations at this 
galactocentric radius can be observed only in the Heavy Sgr 
simulation. A dependance of (Z) with azimuthal angle 6 is 
present. Panel D shows (Z), as a function of R, at 6 — 1.06 
rad. In both discs, a well-defined wave-like pattern can be 
observed, with an amplitude increasing as a function of R. 
Note that (Z) exhibits a stronger dependance with galacto- 
centric radius than with azimuthal angle; a result that could 
be contrasted against currently available stellar samples. In 
the third column of Fig. [4] we explore the relationship be- 
tween these vertical patterns and the radial modes shown in 
Fig.[l] Here, the different colours indicate different values of 
(Z), whereas the relief traces overdense regions. Overdense 
features can be found both above and below the midplane. 



Furthermore, we can appreciate how the mean height of a 
given spiral arm changes as a function of the azimuthal an- 
gle. 

If the observed patterns are indeed signatures of vertical 
density waves, a correlation between (Z) and the mean ver- 
tical velocity (vz) is expected. We explore this in panel E of 
Fig.js] where we plot, for the Light disc, (Z) and (vz) as 
a function of R, at a fixed azimuthal angle. A wave-like pat- 
tern is also observed in (vz). At galactocentric radii where 
(Z) takes an extrema, (vz) ~ 0. On the other hand, at galac- 
tocentric radii where (vz) takes a maximum or a minimum 
value, (Z) ^ 0. This is exactly what is expected from oscil- 
latory behaviour. In the last column of Fig. [4] we show, with 
different colours, the local (vz) for both galactic disc. The 
relief in these panels traces the corresponding overdensity 
maps. Note again the wave-like structure of these patterns. 
In the Heavy Sgr simulation (top panel), (vz) can depart 
from km s""*^ by more than 8 km s""*^. Similarly to what is 
observed for (Z), the (vz) of a given spiral arm varies as a 
function of the azimuthal angle, changing from positive to 
negative departures with respect to (vz) = km s""*^. 



5 DISCUSSION 

In this Letter, we have explored a scenario in which Sgr is 
the perturber behind the North-South asymmetry recently 
observed in the number density and mean vertical velocity 
of Solar Neighbourhood stars (W12). For this purpose, we 
have searched for both local and global signatures of vertical 
density waves in two simulations modelling the response of 
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Figure 5. Panels A and B: Maps of the simulated galactic disc's, 
(Z), at present-day configuration obtained from the Heavy (left) 
and Light (right) Sgr simulations in polar coordinates (in kpc). 
The black dashed lines cross at the location of the blue volume 
shown in Fig^ Panels C and D show the variation of the mean 
height of the disc as we move across the black dashed lines shown 
in the top panels. Blue and red lines are associated with the Heavy 
and Light Sgr simulation, respectively. Panel E: Comparison of 
the mean height (red) and vertical velocity (green), as a function 
of galactocentric radius obtained from the Light Sgr simulation. 



the MW to the infall of Sgr. Distributions of stellar particles 
as a function of height in Solar Neighbourhood-like volumes 
extracted from our more massive Sgr's progenitor simulation 
present clear indications of a perturbation in the vertical di- 
rection of the disc. This asymmetry becomes more evident 
when comparing to a model of the underlying smooth ver- 
tical distribution of particles. Within the |Z|-range allowed 
by our finite mass resolution, the phase-space distribution 
of certain SN-like volumes can qualitatively reproduce the 
North-South asymmetry observed by W12. Remarkably, the 
same phase-space distributions can simultaneously repro- 
duce the signatures of radial density waves observed by G12. 
By creating maps of the mean height of the disc within 



R ^ 20 kpc, we have shown that the vertical perturbations 
observed in local volumes are signatures of a global mode 
perturbing the entire disc. As in the case of the radial den- 
sity modes (G12), the amplitude and the extent to which 
vertical modes can radially penetrate into the disc depends 
on the mass of the perturbing satellite. Interestingly, we have 
shown that the mean height of the disc is expected to vary 
much more rapidly in the radial than in the azimuthal di- 
rection. Furthermore, the mean height of overdense spiral 
features vary azimuthally, moving from below to above the 
midplane of the disc. Signatures of vertical modes should 
also be observable in maps of the Galactic disc's mean ver- 
tical velocity since, not surprisingly, they present a clear 
oscillatory behaviour. 

In contrast to radial modes that can be excited by a 
number of different external and internal mechanisms, ver- 
tical modes must be excited by some agent external to the 
disc. Contrasting the results presented in this work against 
currently available samples of Galactic disc stars could help 
to understand the origin of the observed vertical and in- 
plane perturbations. 
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